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Abstract 

In this work we study the permanence and extinction of a regime-switching predator- 
prey model with Beddington-DeAngelis functional response. The switching process is used 
to describe the random changing of corresponding parameters such as birth and death rates 
of a species in different environments. Our criteria can justify whether a prey die out or not 
when it will die out in some environments and will not in others. Our criteria are rather 
sharp, and they cover the known on-off type results on permanence of predator-prey models 
without switching. Our method relies on the recent study of ergodicity of regime-switching 
diffusion processes. 
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1 Introduction 

In ecosystem, all species evolve and compete to seek resources to sustain their existence. Denote 
the two population sizes at time t by Xt and IJ. Xt denotes the population size of the prey 
and Yt stands for the population size of the predator. The Kolmogorov predator-prey model is 
a general deterministic model taking the form: 

(Xt = Xtf{Xt,Yt), 

\Yt = Ytg{Xt,Yt). 
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For f{x,y) = b — py and g{x,y) = cx — d, one gets the well-known Lotka-Volterra model. 
These deterministic models have been extensively studied. We refer to the monograph due 
to Hofbaner and Sigmund for related study on this deterministic model. Here f{X,Y) and 
g{X, Y) stand for the capita growth rate of each species, which is dependent on the population 
sizes of both species. Various functional response functions have been used considerably in 
modeling population dynamics such as Holling types, Hassel-Varley type, Leslie-Gower Holling 
II type, Beddington-DeAngelis type, etc. However, due to the continuous fluctuation in the 
environment, the birth rates, death rates, carrying capacities, competition coefficients and all 
other parameters involved in this model exhibit random fluctuation to a great extent. To describe 
this phenomenon, stochastic predator-prey models with different kinds of responses are proposed 
and there are many works to study these models (cf. [T], [9], [10] and references therein). For 
a predator-prey system with Beddington-DeAngelis functional response, in practice, we usually 
estimate the birth rate of the prey and death rate of the predator by their average values plus 
errors which follow normal distributions. As a result, one gets the following stochastic predator- 
prey model with Beddington-DeAngelis functional response: 

ciYt 


dXt = Xt{ai-biXt- 


mi+m2Xt+m3Yt 


dt -|- aXtdBi{t), 


dYt = rt[-a,-b,_Y,+ di + PYdBoM. 


( 1 . 1 ) 


with Xq = xq > 0 and Tq = 2/o > 0. All the parameters in (jl.ip are positive, and {B 2 {t)) 

are Brownian motions on the line. For the model (jl.ip . Ji-Jiang [9| and Du et al. [6] studied 
the permanence and ergodicity, Liu-Wang m discussed stochastically asymptotic stability, and 
Li-Zhang [T3| investigated non-persistence and stong/weak persistence. 

From the viewpoint of biological modeling, variability of the environment may have impor¬ 
tant impact on the dynamics of the community. For instance, the distinctive seasonal change 
such as dry and rainy seasons are observed in monsoon forest, and it characterizes the vegetation 
there. Also, in boreal and arctic regions, seasonality exerts a strong influence on the dynamics 
of mammals. Moreover, the growth rates, the death rates and the carrying capacities often vary 
according to the changes in nutrition and food resources. All of these changes usually cannot be 
described by the traditional deterministic or stochastic predator-prey models. Therefore, it is 
natural to consider the predator-prey model in a random environment, which is formulated by an 
additional factor process. More precisely, consider the following regime-switching predator-prey 
model with Beddington-deAngelis functional response: 


'dX, = X,(a,(A,) - (,.(A,)X, - + a{A,)X,dB,(t) 

( - - b,(A,)Y + + l^(Y)y,dB2(t) 


dY = Y 


( 1 . 2 ) 


with Vo = xo > 0 and Yq = yo > 0, where (Bi(t}) and (B 2 (t}) are Brownian motions on the 
line, and (A^) is a continuous time Markov chain with a finite state space S = {1,2,... ,X}, 


2 






1 < < oo. Throughout this paper, the processes {B 2 {t)) and (At) are defined on a 

complete probability space {Q, , {^t}t>o,^), and (At) is independent of {Bi{t)) and {B 2 {t)). 

The parameters ak{-), bk{-), Cfc(-) for k = 1,2, and m/(-) for I = 1,2,3, are all positive functions 
on S. 

The dynamic system (jl.2p is a regime-switching diffusion process, which has been widely 
applied in control problems, storage modeling, neutral activity, biology and mathematical h- 
nance. We refer the readers to [21 n di m [211 [m |23] and the monographs [mils] for the 
study on recurrence, ergodicity, stability, numerical approximation of regime-switching diffusion 
processes with Markovian switching or state-dependent switching in a finite state space or an 
infinite state space. There is a vast literature on population dynamics with regime switching. 
For instance, Du-Du [5] described the omega-limit set of Kolmogorov systems of competitive 
type under the telegraph noise and investigated properties of stationary density; Zhu-Yin [27] 
examined certain long-run-average limits, and Zhu-Yin [28] investigated long-time behavior of 
sample paths for competitive Lotka-Volterra ecosystems. In the study of a population system, 
permanence and extinction are two important and interesting properties, respectively meaning 
that the population system will survive or die out in the future. Yuan et al. [2^ discussed 
extinction for stochastic hybrid delay population dynamics with n interacting species, and Li et 
al. [H] and Liu-Wang [T2| studied permanence and extinction for stochastic logistic populations 
with single species. 

The regime-switching predator-prey model can describe a very important and interesting 
situation. We consider a simple example to introduce it. Let us consider the case S = {1,2}, 
where “1” denotes the rainy season and “2” denotes the dry season. It is rather possible to occur 
the situation that ai(l) — > 0 and ai(2) — ^0:^(2) < 0. This means that the birth rate 

with perturbation of {Xt) in the rainy season makes sure that (Y^) will not die out, but that of 
{Xt) in the dry season makes {Xt) to die out. Then a natural question arises: will (Xt) in model 
o die out or not? This question is very interesting and represents an essential advantage 
of model than model However, the solution of this question is not easy, and so far 

there is few result of this type on the regime-switching predator-prey model or regime-switching 
population dynamics. On this topic, we refer the reader to [TH] for the explicit examples to see 
the complexity of the regime-switching diffusion processes, and to [a da [23 Ea E5] for some 
solutions of this type on ergodicity and stability of regime-switching diffusion processes. In this 
work, we shall provide a sharp criterion to justify whether (Xt) will die out or not (see Theorem 


ll.ll belowL 


Let {qij)ij<=s be the Q-matrix of the process (A*) which means that 



(1.3) 
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for sufficiently small 5 > 0. Throughout this work, the matrix Q = {qij) is assumed to be 
irreducible and conservative, i.e. qtk = —qk ■= — Ylj^k <0. As <S is a finite set and (qij) is 
irreducible, the theory of Markov chains tells us that (A^) is ergodic and there exists a unique 
stationary distribution (fii) for it. To state our main result, we need to introduce two auxiliary 
processes. Let 

dipt = iptiai{At) - bi{At)ipt)dt + a{At)(ptdBi{t), ( 1 . 4 ) 

and 

dipt = a2{At) + - b2{At)-ipt)dt + P{At)'iptdB2it) ( 1 . 5 ) 

V rn2[At) / 

with ipQ = Xq = xq > 0 and tpo = Yq = yo > 0, where (A^), {Bi{t)), {B 2 {t)) are dehned as in 
o. By comparison theorem of SDEs (cf. [8]), ipt > Xt and ipt ^ Yt a.s. for all t > 0. Our 
main result of this work is the following theorem. 

Theorem 1.1 Let {Xt,Yt, At) be defined by Ill.S)) and 11.3)) and {yi)i^s be the stationary dis¬ 
tribution of the process (At). 


{i) If yi(ai{i) — < 0, then lim At = 0 a.s., lim It = 0 a.s.. 

^ ^ y ' t—too t—too 


i£S 


too 


(ii) If y^/ri(ai(i)- cfi{i)) > 0, then {(pt,At) is positive recurrent with stationary distribution 

ies ^ 

TT^ on M+ X S. Assume further 


1 


A - '^yi{a2{i) + -0^{i)) + 


C2l^jX 


i&S 


ies 


/r+ mi(i) + m 2 (i)x 


TT—7r‘^(dx, i) < 0, (1.6) 


then lim Yt = 0 a.s., lim sup At > 0 a.s., and the distribution o/(At,At) converges weakly 

t-too 

to Tr‘^. 

(in) If iai(a 2 {i) H— fi'^ii) -< 0; then {'ipt,At) is positive recurrent with stationary 

V 2 m2(^)/ 

ies ^ ' 

distribution vr’^. Assume further that y^|Uj(ai(i)-ct^(i)) > 0 and 

i€S ^ 


X := X + '^iai(a2{i)+ ^l3‘^{i) -+J2 [ b 2 {i)yTT^{dy,i) > 0. (1.7) 

^ ^ 2 m2ii)J 

r/ien lim sup At > 0 a.s., lim sup It > 0 a.s., and (At,At) At) has a stationary distribution. 

t^OO t^OO 
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This theorem will be proved in next section. As we mentioned, there are few results on per¬ 
manence for the model (|1.2p although there have been numerous works on extinction. Whereas, 
in Theorem 11.11 we provide a criteria which can justify whether a prey die out or not when it 
will die out in some environments and will not die out in other environments. We shall note 
that, when S contains only one state, and hence there is no switching in ()1.2I) in this case, 
our result will coincide with the results in [6]. Actually, according to a similar calculation of 
IS (2.3)], “ 3 ^) + EiG5/R+ & 2 (i)yvr^(dy,i) = 0 in the case A^ = 1, 

hence A = A in this case. But when N > 1, we have no way to calculate the precise value 
of this term at present stage due to lack of explicit representation of the invariant measures 
of regime-switching diffusion processes. Moreover, we should point out that m investigated 
strong (weak) permanence in the mean for the model (II. ip . Compared with [6], the main diffi¬ 
culty in present work is to determine the recurrent property of stochastic processes and 

(V'i, At), which is overcome by using the recent results in [20] on justifying the recurrent prop¬ 
erties of regime-switching diffusion processes. Besides, note that one more condition is needed 
in |6l Theorem 2.2]. Indeed, in the case A > 0, another condition 02 -|- ^ ^ should be 

added to guarantee the process (V’t) in |6] to be ergodic and hence their estimate (2.7) holds, 
i.e. limsupt_^oo j\ny{t) < 0. Without this condition, the estimate (2.7) in [6] is not correct. 

The distributions and are stationary distributions of one-dimensional regime-switching 
diffusion processes. In [2T], for one-dimensional regime-switching diffusion process, an explicit 
representation of the stationary distribution is provided based on the hitting times of this pro¬ 
cess. 

At last, according to the argument of our main result, the permanence or non-permanence 
of (Aj,y)) does not depend on the correlation of the Brownian motions and {B 2 {t)) in 

the situation studied in this work. So in present work, we do not assume any condition on the 
dependence between [Bi{t)) and {B 2 {t)). 


2 Proof of the main result 

We first investigate the properties of the processes and (V't,At) including the estimate 

of their moments and recurrent property. Set 

di := minai(i), ai:=maxai(i) (2-1) 

i£S i£S 

and similarly we can define 0 , 2 , 0 , 2 , 6fc, 6fc, Ck, Ck for k = 1,2 and mi, fhi for I = 1,2,3. By the 
finiteness of S, all the parameters appeared here remain to be positive. 
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Lemma 2.1 For any p > 1 


EipP < (£(/?{]) p exp (^ - (^ai + 


+ 


h 


- I p—L - 2 

ai + 


j-—I 1 - exp I - I ai H- —a 


p — 1 


-p 


( 2 . 2 ) 


and 


EV’f < (E-fAo) ( “ ( “ “2 + ^ + 


+ 


^2 _____ ( ( , C2 , p - 1 -2 


-^2 + + 2^/3^ 


__^l_exp(-(-a2 + ^ + ^/3 )t 


-p 


(2.3) 


Therefore, for any p > 1, 


limsupEpf < (--jzTZJ^ ^ 

t^co ^ai + ‘-^a^' 


and for any p > 1 such that —02 + ^ + > 0; H holds 


(2.4) 


limsupEV’f < 


t—>-oo 


62 


-“2 + ^ + 


-p 


(2.5) 


Proof. We shall only prove the estimate for ipt since the estimate for fjt can be done in the 
same way. By Ito’s formula, 

dpf = p{ai{At) + ^-^a^(Ai)^(pfdt -p6i(Ai)(pf''^Mt +pa(At)(pfdBi(t) 

<p{ai + ^ ^ ^ 6?^iff^dt -pbiipff^^dt+pa{At)(fflABi{t). 

Taking the expectation on both sides and utilizing Holder’s inequality yields that 

<p(ai + - p6iEv9^+^ 

<p{di + ^^-^a‘^^K(pP-pbi{K(pP)~. 

By the comparison theorem of ordinary differential equation, we get (j2.2h , and then (12.4p follows 
immediately. ■ 

By the Lemma [2.11 the processes {(ft, At) and (V’t,At) are nonexplosive, and hence the 
processes {Xt,At) and (1), A^) are also nonexplosive by the comparison theorem. 
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Lemma 2.2 It holds that 


F{Xt > 0, Vt > 0) = 1 and P(yt > 0, Vt > 0) = 1. (2.6) 

Proof. We only consider the process (W) and the result for (Y)) can be proved in the same 
way. Let 

ta = inf{t >0; Xt < A}, A > 0; ax = inf{t >0; Xt > K}, K > 1. 

Let To = inf{t > 0; At = 0}, then ta t as A 0. As (At, At) is nonexplosive, limi^_>.oo = oo 

with probability 1. By Ito’s formula, for p > 0 such that —oi + ^ + > 0, we have 

ftAafcAr^ , -n 4- 1 

+ ( - ai(A^) + 

+ -TWvW + bi{As)Xs) X-Pds 

mi{As) + m2{As)Xs + m3{As)Ys J 

<Xo^ + pE / (-ai + - -h -\-biK)X~^ds. 

Jo ^ rris 2 J 

By Gronwall’s inequality, 

^^ATtXtfATA ^ a^/exp (p( “ «i + ^ + (2-7) 

If P(to < oo) > 0, we can choose t, K large enough so that P(ro < t A aK) > 0. Then for any 
A > 0, 

0 < P(to < t a ax) < P(ta < t a ax) 

T^^KtAafc] — 

< A^XQ^exp - ai + ^ ^ ^ + biK'jt^. 

Taking A | 0, we get a contradiction. Hence P(to < oo) = 0 and we complete the proof. ■ 

Next, we go to study the recurrent properties of the processes (pt, At) and {Jit, At), which 
plays the fundamental role in the proof of our Theorem 11.11 

Lemma 2.3 (Key lemma) (i) When > 0’ ^be process {(ft, At) is 

positive recurrent and has a unique stationary distribution which is a probability mea¬ 
sure on (0, oo) X 5. When /^i(ai(*) — \cx^{i)) < 0, the process {(ft, At) is transient. 

(ii) WhenY^it^s hi{o^ 2 {i) + {i) — m^li) ^ ^ process (V’t, At) is positive recurrent and has 

a unique stationary distribution vr^, which is a probability measure on (0, oo) x S. When 
T.i(zs Piia 2 {iJ + ^/S^(0 - < 0; ibe process (V’t,At) is transient. 
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Proof. (1) We shall use the criterion established in [20] to prove this lemma. As the diffusion 
coefficient of (pt is degenerate at the point x = 0, we use the transform Zt = By Lemma 

12.21 this transform makes sense for all t > 0 a.s.. Then Zt satisfies the following SDE: 

dZt = (ai(Ai) - ^a^(At) - 6i(At)e^*)dt + a{kt)dBi{t). 

The recurrent property of {(ft, At) is equivalent to that of {Zt,At). For each i G S, set := 
(ai(i) — Then the operator £/ defined by 

^f{x,i) = f{-,i){x) + '^ qij{f{x,j) - f{x,i)), f G C2(M x <S) 

is the infinitesimal generator of the process {Zt,At). 

Set Pi = ai{i) — ^a^{i), i G S. We first prove the positive recurrence of {Zt,At). Take 
h{x) = x^ and g{x) = |x|. Then 

= 2(ai(i) — — 6i(i)e^ + for |x| > 0. (2.8) 

Note that linia^^+oo —&i(i)e^ = —oo, lim 2 ;_>._oo —= 0. As > 0, there exist 

e > 0 and tq > 0 such that — s) > 0 and 

2{ai{i) — -C(^{i) — bi{i)e^ + -a‘^{i)x~^)x < —2{Pi — e)g{x), for x > vq, i G S. (2.9) 

2(ai(i) — “ bi{i)e^ + -a'^(i)x~^)x < —2(/3j — e)g{x), for x < —ro, i G S. (2.10) 

Since -2J2i(.s gi{Pi — e) < 0 and lim|^|_j,oo/i(x) = 00 , By [20l Theorem 3.1], there exist a 
Lyapunov function V on R+ x <S with V(x, i) 00 as |x| 00 and a constant tq > 0 such that 


{x, i) < —1, V |x| > ro, i G S. 


Consequently, {Zt,At) and hence {ft, At) are positive recurrent. By [25l Theorem 4.3, pp.ll4], 
{Zt,At) has a unique stationary distribution on R x 5. So the stationary distribution of 
{ft, At) is a probability measure on (0,oo) x S. 

Next, we prove the transience of {Zt, At) under the condition that X]ie<s ^ prove 

it directly due to the behavior of Zt is quit different on (0, + 00 ) and (— 00 ,0) caused by the term 
-6i(z)e*. 

Let h{x) = and g{x) = h'{x) for |x| > 0. Then 


L^'‘^h{x) = (ai(i) — -a^{i) — bi{i)e^ -^)~ 2 ’ for x < 0, i G S. 


X 



As 

lim ai{i) - Iro^ii) - 6i(z)e^ - = j3i, 

x—^—oo I X 

there exist e > 0, ri > 0 so that 

fj,i{/3i + e) < 0, and ai(i) — -a^ii) — 6i(i)e^ — Mx < —ri, i £ S. 

2 X 


As Eie5 + s) <0, the Fredholm alternative (cf. [171 pp.434]) yields that there exist a 
constant k > 0 and a vector ^ such that 


Qi{i) = -K- Pi - e. 

Set V{x,i) = h{x) + iig{x) for x < 0. We have 

s^V(x, i) = LW/i(x) + ^iL^'^">g{x) + Q^{i)g{x) 
+ £ + ii ^ +Q^{i)]g{x) 


^ I Hi 


— ( — K + p 


g{x) 
L^^gjx) 
9{x) 


for X < — ri, i & S. 


It is easy to see that lim 3 ;_»._oo = 0. Denote by ^min = min{^j; z € 5} and ^max = 

max{,^j; i € 5}. Hence, there exists a constant r 2 with r 2 > ri > 0 such that h{x) +Cmin5'(a^) is 
an increasing function on (—oo, —r 2 ] and 


£/V{x,i) < 0, for X < —r 2 , i & S. 


Now, take Zq = z < —r 2 < 0 such that 

h{z)+^ 

max g{z) < h{—r 2 ) + |minfl'(“^2)- 


( 2 . 11 ) 


Set 

T-k = inf{t > i)]Zt = —AT}, r = inf{t > 0; Zt = —r 2 }. 
By Dynkin’s formula, we have 


Ar 

®^[^(-^tAr_xAr ? -^tAr-xAr)] = H(z,io)+E/ s^V{Zs,Ks)<^s 

JO 

< D(z,io), 

where Aq = zq- Letting t —)> oo, we get 

¥.[V{-K,Kr_^)lr_^<r]+W{-X2 1 ^T')^T<r—K ] < V{z,iQ), 
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which yields further that 


{h{-K) + ianng{-K))¥{T > T_k) + {h{-r2) +Cmin5'(-?’2))IF’(T < T-k) < V{z,io), 
p. > N > hjz) + gmaxff(^^) - h{-r2) - gminff(-r2) 

“ ~ H-K) + ^mmg{-K) - h{-r2) - ^mmg{-r2) 

where in the last step we have used the increasing property of x i-)> h{x) + ^ramg{x) on (— 00 , —r 2 ] 
and ()2.1ip . Due to Lemma im —>■ 00 a.s. as K ^ 00 . Consequently, passing to the limit 

as iL — 00 , the previous inequality yields 


P(r = 00 ) > 0, 

which implies that the process {Zt,At) is transient, hence is transient too. 

(2) The results of {'ijjt,At) can be proved by the method analogous to that used above for 
the process (" 0 *, At), which is omitted. ■ 

Proposition 2.4 (Strong ergodicity theorem) Assume that ((/?t,At) and (V’tjAt) are posi¬ 
tive recurrent. Let f be a bounded measurable function on M x 5. Then almost surely 

lim j / f{ips,As)ds = '^ f f{x,i)7r‘^{dx,i), ( 2 . 12 ) 

t Jo Jr+ 

and ^ 

limy/ f{'ips,As)ds = '^ [ f{y,i)TT'^{dy,i)- (2-13) 

t Jo Js.+ 

Moreover, 

liminf- / 6 i(As)<t 9 s ds > [ bi{i)x7r'^{dx,i), a.s., (2-14) 

t Jo Ju+ 

and ^ 

liminfy / 62 (As)V’s ds > / b 2 {i)y {dy,i), a.s.. (2-15) 

t Jo Jr+ 

Proof. We shall only prove (I2.12p and ()2.14p for the process {(ft, At) since (I2.13P and (I2.15P can 

be proved in the same way. The following idea goes back to that of j24l Theorem 3.16, pp.46]. 

For {x,i) G (0, oo)x<S, let 


■= inf{t > 0;((^t,At) = {x,i)}. 

For any {y,j) € (0, oo)x5, (y, j) 7 ^ (x, i), as {(ft, At) is positive recurrent, we have lEj-j, < 00 

and E(yj)r( 3 , j) < 00 . Denote by Ci the time of first return to (x,i) after hitting {y,j), and let 
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(2 = ^CiCi) i-e. C 2 is the time interval between the first time (x,i) is hit after (y,j) has been 
visited and the second such hit. Denote by (n = 0^n-i. Ci, n > 3. By the time-homogeneous 
and strong Markovian properties of {(ft, At), we get that the variables (Cfc) are independent and 
identically distributed stopping times. Moreover, £( 3 , = IE( 3 , -|-E(y < 00 . 

Assume {ifo, Aq) = (x, i) in the following. For a bounded measurable function / on (0, oo)x5, 


/■ELi Cfc ^ f 

/ As)ds = V / 

Jo Tw dy 




= 1 (k 


f{(Ps,As)ds. 


Let 


''Z]fc=i Cfc 


Vi ■= 


J2k^l ^k 


f{iPs,As)ds. 


Then (rji) are mutually independent, and further 


^{x,i)h\ < ll/P(x,i)Ci < 00 , where ||/|| := sup|/(z,fc)|. 

{z,k) 

By virtue of the strong law of large numbers, 

1 n 1 ^ 

lim - Vr/; = E( 3 , j)ryi, lim = E( 3 . i)Ci, a.s. 

n—>-cx) Ti ^ ^ ^ n—>-oo n ^ ^ ' 

1=1 1=1 

Let € N be such that Ci + • • • + Ct < t < (i + ... + Cut+i- Then, for bounded nonnegative 
measurable function /, 


yt yt+i , t't+i vt 

i=l 1=1 /=! Z =1 

It holds that —>■ 00 as t —)• 00 almost surely. Hence, the left-hand and right-hand sides of 

previous inequality both tend to almost surely. Since 

lim j / E,^ i)f{ips,As)ds = '^f f{z,k)TT‘^{dz,k), 

Jo t 7 iJR+ 


it must hold that for any bounded nonnegative measurable function / 


do 


f{z,k)TT‘^{dz,k) a.s.. 


(2.16) 


By the linearity of / in (j2.16p . it is easy to show that (j2.16p holds for all bounded measurable 
function /. So we have proved (I2.12[) . 
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Applying (|2.16p to function f{z,k) = bi{k)'nim{z, M} for positive constant M, we get 
lim - [ 6i(As) min{(/Js, Afjds = [ bi{k)m.m{z, M}7r'^{dz,k) a.s.. 

t Jo 


Then 


1 /■* 1 /■* 
liminf— / 6i(As)(/?5ds > lim — / 5 i(As) min{<y9s, M}ds 
t — ^OO t> Jq t —^OO t, Jq 


keS'^^- 


bi{k) min{z, M}7r‘^(d2;, k) 


a.s.. 


Letting M tend to oo in the previous inequality, we obtain (j2.14p . The proof is completed. ■ 

After the preparation of above results on the auxiliary processes (ipt) and (ipt), we are ready 
to prove our main result. As the proof of Theorem 11.11 is a little long, we divide it into three 
propositions. 


Proposition 2.5 If //i(ai(i) — 
ies 

Proof. By Ito’s formula, we have 


-a^{i)) < 0, then almost surely lim Xt = 

2 t—^OO 


0 and lim 1^ = 0. 

t—>-oo 


dlnX, = („,(A,) - la^(A.) - i„(A,)A, - 

< (ai(At) - ^a^{At))dt + a{At)dBi{t) =: dln<^t. 

By comparison theorem of stochastic differential equation (cf. [8]), we have 

In At < ln(,5t for any t > 0 a.s.. 


Consequently, 


In At ^ in (fit 


In xo 1 /"S /A 

—nXw*- 


^a^(As))ds + y [ a(As)dBi(s). 
^ t Jo 


By the strong ergodicity theorem, 

{t / +7 / a(As)dBi(s)| = y^/it(ai(i) 

-^0 i£S 

Here we have used the fact 



a.s.. 


1 /■* 

lim - / a(As)dHi(s) = 0 
t-tOO t Jq 


a.s., 
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which is due to the boundedness of {cx{k)) and the strong law of large numbers (cf. |15l Theorem 
1.3.4]). Therefore, when < 0, we have limsup^^oo < 0 a.s., which 

further implies that lim = 0 a.s. 

t^OO 

With lim^^oo W = 0 a.s. in hand, we shall show that it must hold lim^^oo Y) = 0 a.s.. 
Namely, when the prey (Xt) dies out, the predator (Y)) must die out too. By Ito’s formula, it 
holds 


, InTt 

lim sup- < 

t—^OO t 


limsup + 1 f f 
t^OO L t t Jq V 


-a2[As) -— 


C 2 (.^s)^s 

mi (As) 


ds + l f l3iAs)dB2is) 
* Jo 


^/ri(o 2 (i) + -/3^(i)) < 0, a.s.. 

ies 


which yields immediately that limt_>.ooht = 0 a.s.. The proof is therefore completed. ■ 

Proposition 2.6 If Hi[ai{i) — -a^(i)) > 0, and A < 0, where A is defined by /il.6\} . then 


ieS 


almost surely limi^oo = 0, limsup^^goW > 0, and the distribution of {Xt, At) converges 
weakly to . 


Proof. We first show that (Y)) tends to be extinct. By Ito’s formula, we have 
1 ..o. C2{At)Xt 


dlnPi = ( - a2{At) - -fi\At) - h2{At)Yt + 


mi(Ai) + m2{At)Xt + mfiAt)Yt 


dt + fi{At)dB2{t). 


Since Xt < g^t a-s. for any t > 0 and {ipt,At) is positive recurrent thanks to Lemma 12.31 by the 
comparison theorem and Proposition 12.41 we obtain 


lim sup — In Yt < 
t—>-00 t 


1 / 1 

limsup - ( - a 2 {As) - -fi^(As) + 

t—loo t Jq \ Z 

1 /■* 

+ limsup - / /3(As)dB2(s) 


C2iAs)<fs 


mi (As) + m 2 (As)v?s 


ds 


t—loo t 

= A < 0, a.s 


Therefore, limi_j.oo Y) = 0 a.s.. 

Next, as limi_).oo A* = 0 a.s., for any e > 0 there exist a measurable subset C 12 with 
P(n£) > 1 — e and a positive constant t(e) such that for any t > t{e), oj E lie. 


dw > ^i(ai(At) - bfiAt)Xt - ^ + ^{^t)XtdBfit) 

V mi(At)+ m2(At)W+ "i3(Ai)e/ 

dXt < Xt{aiiAt) - bi{At)Xt)dt + a{At)XtdBfit), 
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which yields that the distribution of (Xt, A^) converges weakly to the stationary distribution 
of 

To complete the proof, we also need to show that limsup^^ooXi > 0 a.s.. Indeed, if 
^(limj^oo Xt = 0) > 0 , there exists a measurable subset Qq of such that P(no) > 0 and 
Vo; € rio, limt^cxD = 0. For any 5 > 0, as the distribution of {Xt,At) converges weakly to 

vr*^, we have 


7r‘^([0,(5] X 5) > limsupP(Xf € [0,<5]) > P(Qo) > 0. 

t^OO 


By the arbitrariness of 6, we have 7 r‘^({ 0 } x 5) > P(flo) > 0 which contradicts with the fact that 
Tr‘^ is a probability measure on (0, oo) x S (see Lemma [2l^ . We get the desired conclusion. ■ 

Proposition 2.7 If > 0, o,nd \> 0, 

2 2 moil) 

i&S ies ^ ' 

where A is defined by then limsupXt > 0 a.s., limsupYl >0 a.s. and {Xt,Yt, At) has a 

t—)-oo t—>-oo 

stationary distribution. 

Proof. Since almost surely Xt < (ft and Yt < fit for any t > 0, the strong ergodicity theorem 


yields that almost surely 

lim sup — In Xt < lim sup — In 
t—too t , . - e 


t—too t 


1 


“ X] / bi{i)xTT‘^{dx,i), (2.17) 
ies ies 


lim sup — In Yj < lim 

t—too t 


supjln'i/’f + 

t-too t firl ^ m2ii)' 

iGS 

-Yf 

i£S- 


i€S 

-^2 I b2{i)yTr'^{dy,i). 


(2.18) 


Next, we apply the trick used in [ 6 l Theorem 2.2] to estimate the lower bounds. 

1 1 1 /"* 1 

- InXt = - Inxo + - / {ai{As) --a^{As)-bi{As)(fs)ds 
t t ^ Jo ^ 

+ i [\bfiAs){^s - Xs) - iA + 7 /Ws)dBi(s) 

t Jq \ mi{As)+m2{As)Xs+mfiAs)YsI t Jq 

^ilnxo + y/ (ai(A5) - ia^(As) 

t t Jq I 

rt 


bi{As)(ps)ds 

Jo ^ 

t Jo V mi(As) / t Jo 


a.s.. 
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Combining this with (j2.17l) . we get 

liminf j f (bi{As){Xs - (fs) + 
*—>■00 t Jq \ 


ciiAs)Ys 

mi{As) 


^ds > 0, a.s., 


and further 


lim inf — 
t—>oo t 


i(Xs - y?s) + ^ys)ds > 0, a.s.. 


For the process (Yt), we have 

‘ 62 (A.)nds (a 2 {As) + 

C2(As)(/3s C2(A5)^s 


1 1 1 

- In F+ = - In Vo- 

t ^ t ^ t 


C2{As)ips 


1 


t Jo \mi{As)+ m 2 {As)'.ps mi{As) + m 2 {As)Xs 


mi{As) + m2{As)ips 
ds 


(2.19) 


ds 


1 


C2 (As)^s 


C2{^s)Xs 


+ 


t Jo ^miiAs) + m 2 {As)Xs mi{As) + m 2 {As)Xs + m3{As)Ys 
1 


ds 


t 


/3(A,)di?2(s) 


1 


1 


>-\nyo + - - a 2 (As) --/3^(As) + 


C2{As)ips 


C2(As)((^s Xg 


+ 


t Jo ^rni{As) + m2{As)ips \mi{As)m2{As) 


mi{Ag) + m2{Ag)(p, 
C2{As)mo{As) 


ds 


+ 62 (As))lA)ds 


+ 


It follows that 


^ f P{Ag)dB2{s), 
^ Jo 


a.s.. 


lim inf ■ 


1 


C2(As)((/?s Xg 


+ 


C2{Ag)m3{Ag) 

t Jq \mi{Ag) + m2iAg)(pg ' \mi{Ag)m2{Ag) 

‘ ' ■' C2{Ag)(fg 


> 


limsupi f ( -a 2 {Ag)-]-13'^{Ag)- 
t—>oo t Jo ^ ^ 


m\{As)Am2{Ag)kp, 

Invoking ()2.18p and by Proposition 12.41 Lemma 12.31 we get that 


+ 62 (As)jIAjds 

^ds —limsup - Inlt, a.s. 


lim inf — 




+ 


C2{Ag)m3{Ag) 


t^oo t Jq \mi{Ag)+ m 2 {Ag)(pg \mi{Ag)m 2 {Ag) 
> 


+ b2{Ag)^Ygds 


A + ^/ri(a2(i) + + X] hi^Jy{dy, i) = A, 

ies ' 


ies 


a.s. 


and further that 


liminf - / (-^{^fg - Xg) + + fc 2 )h^^)ds > A, a.s 

*->00 t Jq \mi ^ mi 777-2 ' 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 
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Dividing both sides of (|2.19p and (|2.22|] by bi and respectively, and adding them side by 

side, we obtain 


liminfi 

t^oo t Jq C2 Vr?T,2 


+ 


hrhi , Cl \-i- 


+ 


A, a.s.. 


C2 rhibi 

Hence, if A > 0, then inequality (12.231) implies that it must hold limsupi^oo Yt > 0 a.s.. 
Note the following facts: 

C2(At)Xi 


(2.23) 


dlnYt= (^-a2{At)-^/3\At) + 


mi{At) + m2{At)Xt + m‘j,{At)Yt 


dt + /3(At)dH2(t), 


and limt^oo \ Jq /3(As)di?2(s) = 0 a.s.. If P({w : limt_>.oo Xt{io) = 0}) > 0, then it must hold 


P(lim sup - In < o) > 0. 
^ t—^oa t / 


Hence, the fact limsupj_j,oo 1* > 0 a.s. implies that limsup^_^oo > 0 a.s.. 

At last, we show the existence of invariant probability measure of the process {Xt,Yt, At). 
Note that almost surely Xt < (pt and Yt < il^t for any t > 0. By (12.4p and p2.5p of Lemma 12.11 
we have that for p > 1 satisfying —02 + ^ + > 0, there exists a constant C > 0 such that 

E(XP + YP + Af)ds < C. 

According to O Theorem 4.14], there exists a stationary distribution for (Xt, Yt, At). The proof 
is completed. ■ 

Proof of Theorem 11.11 

The argument follows immediately from the arguments of Propositions 12.3112.5112.71 ■ 
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